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We use the Bethe ansatz equations to calculate the charge stiffness D c = 
(L /2)d 2 Eq/ d&c\$ c= Q of the one-dimensional repulsive-interaction Hubbard model for 
electron densities close to the Mott insulating value of one electron per site (n = 1), 
where Eq is the ground state energy, L is the circumference of the system (assumed 
to have periodic boundary conditions), and (hc/e)<f> c is the magnetic flux enclosed. 
We obtain an exact result for the asymptotic form of D C (L) as L — ► oo at n = 1, 
which defines and yields an analytic expression for the correlation length £ in the 
Mott insulating phase of the model as a function of the on-site repulsion U. In the 
vicinity of the zero temperature critical point U = 0, n = 1, we show that the charge 
stiffness has the hyperscaling form D c (n, L, U) = Y+(£<5, £/L), where 6 = |1 — n| and 
Y + is a universal scaling function which we calculate. The physical significance of 
£ in the metallic phase of the model is that it defines the characteristic size of the 
charge-carrying solitons, or holons. We construct an explicit mapping for arbitrary 
U and £<5 <C 1 of the holons onto weakly interacting spinless fermions, and use this 
mapping to obtain an asymptotically exact expression for the low temperature ther- 
mopower near the metal-insulator transition, which is a generalization to arbitrary 
U of a result previously obtained using a weak-coupling approximation, and implies 
hole-like transport for < 1 — n <C 
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I. INTRODUCTION 



The Mott-Hubbard metal-insulator transition M is one of the fundamental problems in 
electronic condensed matter physics. The assertion by Anderson [[| that the high-T c CuO"2 
superconductors should be viewed as lightly doped Mott insulators has stimulated much 
work on the problem. In this paper, we use the Bethe ansatz equations |3H| to calculate the 
zero temperature conductivity and low temperature thermopower of the one-dimensional 
(ID) repulsive-interaction Hubbard model for electron densities close to the Mott insulating 
value of one electron per site. 

One way to characterize interacting fermion systems is via the zero temperature fre- 
quency dependent conductivity a(uj), which is the linear response of the system to a spatially 
uniform, time dependent electric field. In general, cr(u) has the form 

a{u)=2D c (^-+TT8{uj)\+a reg {uj), (1.1) 

where lim ua reg (uj) = 0. The coefficient D c is the charge stiffness. In a Galilean-invariant 

system with n electrons per unit volume of charge e and mass m interacting via velocity- 
independent forces, D c = ne 2 /2m and a reg {uj) = 0, independent of the interactions. For 
electrons moving in an external potential, D c may depend on the interactions and on the 
potential. Walter Kohn has argued || that for spatially infinite systems with a discrete 
translational invariance two cases are possible: either D c = and lim Re[a r . e9 (c<j)] = 0, or 

D c > 0. If D c = 0, the material is insulating; if D c > 0, the material has the infinite 
dc conductivity and perfect diamagnetism expected of a perfect metal. D c is thus an or- 
der parameter |J for the Mott-Hubbard metal-insulator transition which may occur in a 
nondisordered system as the electron concentration and interaction strength are varied. 

In a system of finite length L with periodic boundary conditions, D c is in general non- 
vanishing 0. In a previous paper j7|, we showed numerically and via a large-?/ approxima- 
tion that in the Mott insulating phase of the ID Hubbard model, which occurs at n = 1 for 
repulsive interactions, the large- L behavior of D c is 

D C (L) ~L 1 / 2 exp(-L/0- (1.2) 

This equation defines the correlation length £ of the Mott insulator, which is a function of 
the on-site repulsion U. In this paper, we derive Eq. (|1.2|) analytically for arbitrary U, and 



show that the correlation length so defined is important also in the metallic phase of the 
model. In particular, we show that in the vicinity of the zero temperature critical point 
U — 0, n = 1, the charge stiffness assumes the hyperscaling form: 

D c (n,L,U)=Y + (£5,£/L), (1.3) 

where 5 = |1 — n| and Y + is a (presumably universal ||) scaling function which we calculate. 

The physical significance of £ in the metallic phase of the model is that it defines the 
characteristic size of the charge-carrying solitons, or holons. We show that Woynarovich's 
reformulation M of the Bethe ansatz equations of the ID Hubbard model in terms of the 
parameters of the charge excitations only is formally equivalent to an asymptotic (large 
L) Bethe ansatz for holons, and use this asymptotic Bethe ansatz to construct an explicit 
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mapping for arbitrary U and £5 <C 1 of the low- lying charge degrees of freedom of the model 
onto weakly interacting spinless fermions. Using this mapping, we obtain asymptotically 
exact expressions for the charge stiffness and low temperature thermopower near the metal- 
insulator transition. Our result for the thermopower extends and makes more rigorous 
previous work of Schulz [HJ which was based on a weak-coupling approximation, and implies 
that the transport is hole-like for < 1 — n <C 

The paper is organized as follows: In Sec. O, we review the linear response arguments 
and Bethe ansatz equations used to calculate D c in ID. In Sec. |TTTj, we present an exact result 
for the charge stiffness at n = 1 in the large- L limit, verifying Eq. ( |1.2| ), and show that the 
correlation length so defined also governs the exponential decay of the equal-time single- 
particle Green's function at n = 1 in both the weak- and strong-coupling limits, as well as 
the pairing correlations in the ground state of the attractive ID Hubbard model at n = 1. In 
Sec. |V|, we determine the scaling behavior of D c near the metal- insulator critical point. The 
total optical spectral weight 7riV tot = f£° Re [a{ui)]dui is also discussed. In Sec. [V], we calculate 
the finite-size corrections to D c and N tot in the regime L > £, and discuss the implications 
for the interpretation of numerical calculations of ct(uj) in interacting fermion systems on 
small clusters. In Sec. [VI|, we determine the charge excitations, charge stiffness, and low 
temperature thermopower near the metal-insulator transition. We summarize our results in 
Sec. |V1J| , and comment briefly on their relevance for the physics of high-T c superconducting 
materials. The details of the calculations of D c at and near n — 1 in the large-L limit are 
given in Appendixes A and B, respectively. 



II. FORMALISM 

We consider the Hubbard model of spin-1/2 fermions hopping with matrix element t 
between nearest-neighbor sites of a one-dimensional lattice with unit lattice constant, and 
subject to a repulsive interaction U when two fermions (of opposite spin) occupy the same 
lattice site. We impose periodic boundary conditions, and thread the system with a (dimen- 
sionless) spin-dependent flux which we represent by a uniform spin- dependent vector 
potential A a = (he/ 'e)§ a / L, where L is the circumference of the system. The Hamiltonian 
may be written 

H = ~t £ [e^cj^cy + H.c] +Uj2n n n lb (2.1) 

l,cr I 

where n icr = c\ a Ci a and the sums run from I = 1 to L and cr =|, j. In the following, we set 
e = h = t = 1 unless explicitly stated. 

The numbers M and M' of up- and down-spin electrons are constants of the motion; we 
denote the minimum energy eigenvalue in a given sector (M, M') by E(M, M', U, $f, $|). 
The sign of the interaction U in Eq. (|2.1|) can be reversed by a particle-hole transformation 
on the down-spin electrons [ Tt \ , 



c \i ~^ ( lV c U> c z\ > c lt> (2 2) 



leading to the following relation, 
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E(M, M', U, $ T , $J = E(M, L - M', -U, $ T , + M?7, 



(2.3) 



which we will exploit in Sees. [TTI| and |V1J| . 

To simplify matters, we consider only systems where the total number of electrons N 



M + M' is even. Then the ground state is a singlet [15] and the ground state energy is 
E = E(N/2, N/2, U, $ T , The charge stiffness D c is defined by 



ld 2 (E /L) 



2d($ c /L)< 



(2.4) 



where $ c = (<&f + <&J/2. Applying second-order perturbation theory in $ c to Eq. fl2.1|) gives 



an 



(2.5) 



where T = — X^( c J+1ct c ^ + #-c.) is the kinetic energy operator, J p = —iJ2( c l+i a c iv ~ H.c.) 
is the paramagnetic current operator, ( ) denotes the expectation value in the ground state, 
and all quantities are evaluated at $^ = $| = 0. 

We next consider threading the Hubbard ring with a time-dependent flux $f(t) = = 
$ c (0) exp(— iut), which leads to a uniform time- dependent electric field by Faraday's law. 
Standard linear-response arguments applied to Eq. fl2.1|) give the frequency-dependent con- 
ductivity at zero temperature 



(TlUJ 



(-T) 2 \(v\J P \0)\ 2 
L L^uj-E. + Eo + iO- 



(2.6) 



The conductivity thus has the form (1.1) with D c given by Eq. (|2.5f) . The high-frequency be- 
havior of Im [<t(u;)] leads, via the Kramers-Kronig relations which link the real and imaginary 
parts of cr(u)), to an /-sum rule for the total optical spectral weight @,14], 



TtN 



tot 



J o °°Re[a(u J )}dco=^ L {-T). 



(2.7) 



The Hellman-Feynman theorem gives (T) = E — UdE /dU. Thus both N tot and D c may 
be obtained as derivatives of the ground state energy. 

The energy and momentum of the eigenstates of Eq. (j2.1|) can be expressed via a gener- 
alization of Bethe's ansatz as 



N 



E = — 2 ^2 cos K. 



n=l 



N 



(2.8) 



(2.9) 



71=1 



where the pseudomomenta k n are a set of distinct numbers determined by 
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Lk n = 2nl n + ^ - £ 2 tan" 1 (2 .l ) 
and the spm rapidities A a are a set of distinct numbers related to the k n by || 

E 2 tan- 1 Aa ; S f fc - = 2^ + $ T - ^ + £ 2 tan^ 1 (2.11) 

n=l U / 4 /3=1 U / Z 

Here {J n ; n = 1, . . . , N} and { J Q ; a = 1, . . . , M} are the quantum numbers which specify 
the state of the system |T5|: the I n describe the charge degrees of freedom []IB|,[T7]|| , and 



are distinct integers (half- integers) for M even (odd) ||, while the J a describe the spin 
degrees of freedom |T6|JT7| ,|9[| , and are integers (half-integers) for iV — M odd (even) ||. For 
the ground state, the quantum numbers are consecutive integers (or half-integers) centered 
about the origin ||. 

In the following, we solve the Bethe ansatz equations ( |2. 10| ) and ( |2.11| ) analytically 
and numerically in various limits to obtain D c and N tot as a function of n, L, and U . 
The knowledge of the eigenstates of H also allows us to calculate the low temperature 
thermopower near half filling. 



III. THE CORRELATION LENGTH OF THE MOTT INSULATOR 

A. Definition of £ 

The ID Hubbard model for n = 1 and U > is known to have an insulating ground 
state , so that the charge stiffness D c vanishes for an infinite system. In a previous paper 
0], we showed numerically and via a large-f/ approximation that for a large finite ring of 
circumference L at n — 1 

D C (L) = (-1)* +l L l ' 2 D(U) exp[-L/£([/)]; L - oo (3.1) 

(L even), where D(U) is a positive [/-dependent number. Eq. ( ^.1[ ) serves to define the 
correlation length £,(U) in the Mott insulating phase of the ID Hubbard model. We have 
subsequently derived Eq. (|3.1|) analytically for arbitrary U. The details are given in Ap- 
pendix A. The result for the correlation length £(£/) is 



l/m = uk dy cosh(2ny/U) " (3 ' 2) 



The result for D(U) is given in Eq. ( |A28|) . The limiting behavior of D(U) is: 



^{^ ),/2 S=£ (3 - 3) 

where A ~ 0.147376. 

The negative sign of D c when L is a multiple of 4 indicates that the persistent current 
of the ring is paramagnetic, as discussed in Refs. [0,0 • Orbital paramagnetism is a generic 
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feature of half-filled single-band An electron systems, and was first observed many years ago 
i~9fl in NMR spectra of [16]annulene and [24]annulene, larger analogs of benzene. The full 



flux dependence of the ground state energy of the Mott insulator is 

£ ($ c ) - Eq{0) = - cos $ c ); L — > oo. (3.4) 

By contrast, in the metallic phase at n 7^ 1, the flux dependence of the ground state energy 
in the large-L limit is E (§ c ) - E (0) = D C $ 2 C /L for |$ c | < 2irL, and the periodicity 



Eq(&c) = E (§ c + 2tt) required by gauge invariance p0| is restored by level crossings. 



B. Relation of £ to the Single- Particle Green's Function 



The zero temperature equal-time single-particle Green's function is defined as 

G aa ,(\x -A) = (0|4v<W + cLc,v|0). (3.5) 

From Eq. ( |2.1| ) it follows that G aa i{\x — x'\) = G aa (\x — x' \)S a>(T '. We wish to show that the 
Green's function of the Mott insulator has the form (at L = oo) 



G a(T (\x\) ~ exp(-|x|/f) as 



oo. 



(3.6) 



In order to do so, it is useful first to consider some limiting cases. 

The weak-coupling limit of £ may be obtained from an asymptotic expansion of Eq. ( |3.2| 
around U = 0. We obtain 



lime = ^(t/U) 1/2 e 27Tt/u (l + U/lGnt + 



2t + U/2n + • 



A(U,t) 

where we have made the hopping matrix element t explicit, and 



(3.7) 



2A(U,t) 



16t 2 



(y 2 - l) 1/2 dy 
U J i sinh(27rfy/[/) 



(3.8) 



is the Lieb-Wu charge gap p|,P,[T6|| . The low energy (continuum limit) physics of the ID 
Hubbard model has been studied extensively by bosonization techniques pl|j22]j23|j2^j25[] . In 
the low energy sector, the effective Hamiltonian separates into independent terms describing 
the charge and spin degrees of freedom [ 2lL|22 ,|23 , 2i , 25] . Near n — 1 and in the weak-coupling 



limit, the Hamiltonian for the charge degrees of freedom is of the relativistic sine-Gordon 



form The quantity which plays the role of the speed of light c in the Hamiltonian is 

given to leading order in U at n — 1 by [pl|,25 



c = 2t + U/27T + 



(3.9) 



The elementary objects in the sine-Gordon model are solitons [2S] carrying charge e and 
obeying the relativistic dispersion relation E^ = (c 2 k 2 + A 2 ) 1 / 2 [27,28]. The quantum soliton 



length £ s = c/A defines the characteristic size of the solitons P9"| , so we see from Eqs. ( \$-7\ ) 
and ( |3.9| ) that in the weak-coupling limit the correlation length is equal to the quantum 
soliton length. 
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The mapping onto the weak-coupling Hubbard model fixes the number density of solitons 
as n s = 1 — n [|5[]. The operator c xa involves creating at least one soliton (plus spin 
excitations). Thus, if we evaluate Eq. ( |3.5| ) at a time-like separation at n = 1, we find 
G aa (t—t') ~ exp[— iA(t— t')} + terms involving higher energies. Using the Lorentz invariance 
of the weak-coupling field theory to rotate back to a space-like separation, we obtain G aa (\x— 
x r \) ~ exp(— A|x — x'\/c) [j30| . The gapless spin degrees of freedom correct this expression 
only by a power of \x — x'\. We have thus established Eq. ( ^.6| ) in the weak-coupling limit. 

Now let us consider the strong-coupling limit. From Eq. ( |3.2| ) it follows |31|] that 



r 1 = HU/at); U^oo, 



(3.10) 



where a = (r(l/4)/\/27r) 4 ~ 4.377 (this result was obtained previously by us in Ref. 0). 
Eqs. and O then im P!y that D c(L) = AUL l l 2 {at/U) L , where A ~ 0.147376. The 
correlation length may thus be obtained in the large-?/ limit from an expansion of the ground 
state energy in powers of t/U. The leading contribution to D C (L) comes from the lowest 
order term which is flux dependent, and is M 



lim 

U->oo 



2D C (L) 



t l 



e (0.14 n ( d 



L-l 



where v\ = Eac\ +la c la , 



D 



Pes L 



-i„,t 



|0oo) + E.c. 



(3.11) 



l 'HX'Hh 



|0oo) 



lim^oo |0), and Sl is the group 
of permutations of L objects. This expression for D c is a sum over all processes which 
transport one unit of charge around the ring in the minimum number of steps, L. In the 
large-?/ limit, the Green's function G aa {L — 1) may be evaluated by expanding the ground 
state |0) in Eq. Q to 0((t/U) L ^). The resulting expression for the leading contribution 
to i[G TT (L - 1) + G U (L - 1)] is identical to that for 2D C (L)/L in Eq. (|37TT| ) , except that 
|0oo) now refers to the infinite system, not the periodic system of length L ||32|| . One can fix 
the phase factors such that = G^. Then the only difference between tG aa (L — 1) and 
D C (L)/L in the large-?/ limit stems from the difference between the spin correlations in the 
infinite Heisenberg model and those in the periodic one of size L, which can at most cause 
tG aa [L — 1) to differ from D C (L)/L by an algebraic factor in L. The correlation lengths 
defined by the two functions must be the same in the large-?/ limit. 

We have thus established Eq. ( |3.6| ) in both the weak- and strong-coupling limits. We 
now give a physical argument which suggests that G has this asymptotic form for arbitrary 
U. 



C. Singlet Pairing for U < 

The modulus |G ao .(|a;|)| of the equal-time single-particle Green's function is invariant 
under U — > — U at half filling, as may be seen by applying the transformation ( |2.2|) to 
Eq. (|3.5|) . Let us therefore consider Eq. ( |3.6|) for the case U < 0. The U < Hubbard model 
is a superconductor (or more precisely, in one dimension, has divergent superconducting 
fluctuations p3| ). Superconductors are characterized by a correlation length £ which one 
may think of as the size of a Cooper pair, and G decays exponentially with characteristic 
length £. The superconducting correlation length can also be determined via the periodicity 
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of the ground state energy of a superconducting ring with respect to the magnetic flux 
enclosed. Gauge invariance implies that E must be periodic in $ c with period 2ti (one flux 
quantum) |^D) ; however, in a system with superconducting correlations, the period will be 7r 



(half a flux quantum) when the system is large compared to the correlation length [p0[| . In 



this subsection, we show that the deviations from a periodicity of half a flux quantum in the 
ground state energy of the U < ID Hubbard model scale as exp[— L/£(|Z7|)] at n — 1. We 
therefore identify £(|t/|) as the superconducting correlation length; this makes it plausible 
that G decays exponentially with this length for arbitrary U, as asserted in Eq. ( |3.6| ). The 
analogous electron-hole correlations in the ground state of the Mott insulator were discussed 
by Krishnamurthy et al. in a slightly different context. 

The Bethe ansatz description of the ground state of the attractive ID Hubbard model 
involves complex pseudomometa ||35|| , and is unwieldy for a finite system with periodic 
boundary conditions. However, we can use the identity (|2.3|) to obtain the ground state 
energy of the half-filled attractive ID Hubbard model as a function of $ c from the ground 
state energy of the half-filled repulsive ID Hubbard model as a function of $ s = ($j — $;)/2. 
The latter may be obtained from a solution of Eqs. ( |2.10| ) and ( [2. lip with the fc's and 



A's strictly real (alternatively, one can think of this as a Bethe ansatz for U < in a 
representation where the down-spin electrons are treated as holes). In ID, the variation of 
the ground state energy with respect to flux for U < is E (§ c ) ~ E (0) + D c (— \U\)$ 2 C /L 
when L is large. However, there are level crossings with states whose energy parabolas have 
local minima at $ c = ±7r, ±27r, etc. Each parabola corresponds to a particular choice of 
the Bethe ansatz quantum numbers {/„} and {J a } (in the repulsive representation). The 
quantum numbers {/°} and of the ground state at $ c = are given in Eqs. ( |A1| ) 

and (|A2|) in Appendix A. The ground state at $ c = 2tt has quantum numbers satisfying 
I n = I® + 1, J a = J° — 2, and is equivalent to the above state, but shifted by one flux 
quantum. The equivalence of these two states follows from gauge invariance |p0| . There is 



a level crossing with a third, inequivalent state whose energy minimum lies at <3> c = 7T, and 
whose quantum numbers satisfy /„ = 1°, J a = J° — 1. The energies of these three states as a 
function of $ c are plotted as solid curves in Fig. [I] for U = —4 and N = L = 6 and 28 (note 
that £(4) w 4.06). The ground state energy for the repulsive case (shifted by — N\U\/2) is 
plotted as a dashed line for comparison. Equations Q2.10|) , (|2.11|) , and (|3.4j ) imply that in 



the large-L limit the local energy minimum at $ c = it differs from that at <3> c = by 

E (n) - E Q (0) = „ L-^eM-L/mi)}- (3-12) 

It is evident from Fig. [I] that as L increases and the local energy minima at $ c = and 
$ c = 7r become degenerate, the ground state energy becomes periodic with period 7r (half a 
flux quantum). 



IV. SCALING BEHAVIOR OF D c NEAR THE 
METAL-INSULATOR CRITICAL POINT 

We have shown that the charge response of the half-filled (n = 1) U > ID Hubbard 
model may be characterized by a length £ which diverges as U — > 0. In this section, we 
show that the point U = 0,n = lisa conventional quantum critical point in the sense that 
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the singular behavior of D c (n, L, U) in the vicinity of this point is given by a (presumably 
universal ||) scaling function Y(£|l— n\, £/L), which we calculate numerically and, in various 
limits, analytically. Our results confirm the applicability of the hyperscaling ansatz to this 
system. We also discuss the total optical spectral weight 7rN tot . 

The hyperscaling ansatz for a T = phase transition is that the singular part of the 
ground state energy per correlation volume scales as the inverse of the correlation time || . 
Since the correlation time scales as £ z [ |36|| , it follows that 

E° ing 7 'L d ~ C (d+Z \ (4-1) 

where d is the spatial dimension and z is the dynamic critical exponent. In our problem 
d = 1, and the Lorentz invariance of the critical point implies z = 1. From Eq. (2~%), we 



see that D c involves two derivatives of Eq/ L with respect to & C /L; the singular part of D c 
should therefore have scaling dimension 2 — d — z = 0. We thus postulate that as U — > 0, 
n — *■ 1, and L — > oo, the singular part of D c is given by 

Df»(n,L,l7)=y±(^e/i), (4-2) 

where 8 = \l — n\ and Y± are dimensionless scaling functions, universal up to a metric factor 
fixing the units of length and energy. The subscripts ± on Y refer to U > and U < 0. 
The charge response at £7 = 0~ is that of noninteracting electrons with nearest neighbor 
hopping, so y_(0, 0) = 2/ir. One consequence of Eq. ( |4.2|) , therefore, is that the charge 
stiffness at n — 1 and L = oo has a universal jump of 2/ir as C/ goes from + to 0~. 
Eqs. (O) and (O) imply that 



Y + (0,y) = (2/ny) 1 / 2 exp(-l/y); y ^ 0. (4.3) 
To verify the scaling law (O) for the general case, we have calculated D c (n, L, U > 0) by 



solving the Bethe ansatz equations ( 2 . 1 0|) and ( 2.11|) numerically for systems with < 200 



and L > N. We have considered only systems with an even number of electrons N. In order 
to obtain for D c a smooth function of n and L, we impose antiperiodic boundary conditions 
(i.e., $^ = $j = 7r) when iV mod 4 = 0, so that the ground state is always a singlet |3"T|,|3"g 



Because 5 = |1 — N/L\ = 0, 1/L, 2/L, . . ., the scaling function Y(x, y) is only defined on a 
countable set of lines in the x, y plane, namely on x = my; m = 0, 1, 2, . . . , oo. In Fig. |2|, 
we plot ttD c versus £/L along the lines 5 = 0, 1/L, 2/L, and 4/L for systems with N = 60, 
80, and 100 electrons. In each case, the results for iV = 60, 80, and 100 fall onto a smooth 
curve over the entire range 0.3 < £/L < 10000, verifying Eq. (|4.2|) . (The limit L 3> £ is 
examined more fully in Sec. 0.) For comparison, Fig. 0(b) is replotted with £ as the abscissa 
in Fig. 0(e). It is evident from Fig. |2|that as£/L— >oo, D c ^ 2/tt (the value of the universal 
jump at the critical point), which we may express as 

2 

Jim Y+(x,y) = -. (4.4) 



The charge stiffness in the limit L — > oo (y — >• 0) has been calculated previously [24,39 
from a numerical solution of the Bethe ansatz integral equations. Here we merely point out 
that because Y is universal, one can extract Y + (x,0) from Haldane's analytical result |29| 



for the sine-Gordon model at j3 = 8n . The mapping between the two models is uniquely 
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determined by equating the energy gap A to the soliton rest energy m s c 2 , equating the 
correlation length £ to the quantum soliton length, and equating the doping 5 = |1 — n\ to 
the mean number of solitons per unit length, as discussed in Sec. [T I [T3[ The result for the 
scaling function is |29| 

Y + (x, 0) = x (l - \{kx) 2 + ^(vrx) 3 + ■■); x - 0, (4.5) 



Y + (x,0) = -(l--[\n{nx)]- 1 + ■■■); x -> oo. (4.6) 

Near the critical point, therefore, the charge stiffness grows rapidly as one dopes the Mott 
insulator, saturating near its maximum value 2 /it when x — £5 ~ 1. In Sec. VI, we com- 
pute the small- x behavior of the scaling function directly from the Bethe ansatz equations, 
verifying Eq. ft4.5|) . 

The total optical spectral weight 7rN to t involves information at high frequencies, but is 
related via the /-sum rule ( |2.7|) to the expectation value of the kinetic energy in the ground 
state; thus iV tot may also be expected to exhibit universal behavior near the critical point. 
In the critical regime £, L, ^> 1, we find that pO 



iVtot * -, (4.7) 

7T 

approximately independent of 5, L, and £. The approximate 5-independence of N tot near the 
critical point is in agreement with the intuitive picture that when £ ^> 1 the kinetic energy of 
the system should be little affected by the transition to an insulating state. At 5 = 0, there is 
an energy gap 2A for charge excitations, and all the optical spectral weight is at frequencies 
u > 2 A when L = oo. We refer to the states at ui > 2 A as being in the "upper Hubbard 
band" (UHB) (see Sec. |VJ). Although for 5 > there is no true gap in criuj) = Re [5"(u;)], 
the structure in the UHB region u > 2A must persist for sufficiently small 5 by continuity. 
However, the rapid growth of D c with 5 and the approximate (^-independence of N tot imply 
that near the critical point the primary effect of doping on ct(uj) is to transfer spectral weight 
from the UHB region into the Drude peak at uj = 0, so that the UHB structure in ct(oj) is 
essentially destroyed for 5 > A similar shift of spectral weight from the UHB region 
into the Drude peak with decreasing system size is implied by Fig. the UHB structure 
being essentially destroyed for L < £. In the language of the continuum field theory [^T 



the approximate equality N tot ~ D c ~ 2/7T, which holds in the weak-coupling limit both for 
£5 3> 1 and for £ ^> L, arises because the effects of umklapp scattering are largely smoothed 
out in those limits, so that d-{uj) approaches the result expected in the absence of umklapp 
scattering with N tot = D c = (2/w) sin(n7r/2). 



V. FINITE-SIZE CORRECTIONS TO THE OPTICAL SPECTRAL WEIGHTS 

In this section, we calculate the finite-size corrections (for L > £) to the charge stiffness 
D c and the total optical spectral weight vrA tot in the metallic phase of the U > ID Hubbard 
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model. Since the scaling function for the charge stiffness is expected to be universal || , the 
finite-size corrections we obtain should be relevant for other, nonsoluble models as well. 

We find that the finite-size corrections to both D c and N tot are positive and depend 
smoothly on L, provided the boundary conditions are chosen so that the ground state with 
an even number of electrons is a singlet: periodic boundary conditions when N mod 4 = 2, 



anti-periodic when N mod 4 = |42]. We employ these boundary conditions throughout 
this section. In Fig. |3], D c and N tot are displayed as a function of L for 5 = 0.2 and U = 6 
(£ ~ 2.04). Note that the finite-size corrections to D c are larger than those to N tot ; this is a 
generic feature of the finite-size corrections, and indicates that the finite-size enhancement of 
D c is in part due to a transfer of spectral weight from higher frequencies, as discussed above. 
We find that the leading finite-size corrections to both D c and Ntot are 0(L~ 2 ) in the metallic 



phase |4"5fl, in accordance with the system-size dependence of the ground state energy [04 



The solid curves in Fig. |3| are extrapolations of the L~ 2 behavior to small L, indicating that 
the asymptotic form of the finite-size corrections gives a reasonable approximation even for 
fairly small systems, provided L > £. 

The coefficients A(S,U) and B(5,U) of the leading finite-size corrections in the metal- 
lic phase, defined by D C (L) = D c (oo)[l + A(5,U)/L 2 + • ■ ■] and N tot (L) = iV tot (oo)[i + 
B(5,U)/ L 2 + ■ • ■] are plotted versus U for several values of 5 in Fig. 0j (note the different 
vertical scales). Both for U — *■ and far from the critical point, A(5, U), B(S, U) — ► 7r 2 /6, 
the value expected for noninteracting fermions (with or without spin) with nearest-neighbor 
hopping. While the relative finite-size corrections to the total optical spectral weight are 
always of order tt 2 /6L 2 , the relative finite-size corrections to D c can be much larger, par- 
ticularly near the critical point. The curves in Fig. HJ(a) display increasingly sharp maxima 
as 5 is decreased, which occur at values of U such that £(U)5 ~ 0.4, and at which point 
A{5,U) « 0.8«£(£/) 2 . 

In the critical regime, we expect the leading finite-size corrections to D c to scale as 

D C (L) = D c (oo)[l + n 2 /6L 2 + / + (£6)(£/ L) 2 + •••], (5.1) 

where f+(x) = d 2 Y + (x, y)/dy 2 \ y=Q /2Y + (x, 0), and we have included the regular term 7r 2 /6L 2 . 
The small-x behavior of f+(x) can be obtained by noting that to 0(!;5) 3 , the charge stiffness 
in the weak-coupling limit is equivalent to that of noninteracting relativistic spinless fermions 
H [c.f. Eq. (p])], whence 



/ + (a;) = y-7rV + ...; x - 0. (5.2) 

In Fig. H the function / + = [A(8, U) — n 2 /6}/C,(U) 2 is plotted versus the solid curve 

was obtained by fixing 5 and varying U, while the triangles were obtained by fixing U and 
varying 5. Note that the numerical data connect smoothly to Eq. ( |5.2|) (dotted curve) as 
£5 — > 0. Fig. |^ provides further proof of the scaling law ( [4.2| ) in the regime L ^> ^. 

Though the finite-size corrections to D c near the critical point are much larger than 
expected in a noninteracting system, they should not be a great obstacle to obtaining infor- 
mation about optical properties of ID models from numerical calculations on small lattices, 
except at 5 = 0, due to the impossibility of simultaneously having small L and small but 
finite 5 (min<5 = 1/L). Roughly speaking, the maximum value of the L~ 2 correction in the 
critical regime is 
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D c (L)-D c (oo) _ 0.13 
D c (oo) ~ 

though of course we've seen in Fig. |3| that the additional term 7r 2 /6L 2 can lead to slightly 
larger corrections (~ 20%) for small systems. 

A further comment on the choice of boundary conditions is appropriate at this point. 
One practice which has been employed to reduce finite-size effects is to average over periodic 
and antiperiodic boundary conditions. While this may result in cancellations for some values 
of n and U due to the fact that the finite-size corrections sometimes have different signs, 
it is ineffective near the critical point, where the finite-size corrections are largest, since 
the corrections are always positive and essentially independent of boundary conditions in 
that regime |j45| |. Furthermore, depending on U and n, any of the three possibilities — either 



periodic or antiperiodic boundary conditions or an average over the two — can result in the 
smallest finite-size corrections. One is better off choosing boundary conditions such that the 
ground state with an even number of electrons is a singlet; then the finite-size corrections 
to D c and N to t are always positive, with magnitudes as given above. 



VI. CHARGE EXCITATIONS AND THERMOPOWER NEAR 
THE METAL-INSULATOR TRANSITION 

The mapping of the charge sector of the weak-coupling ID Hubbard model near n = 1 
onto the sine-Gordon model suggests that the charge excitations near half filling can be 
described in terms of the states of L — N soliton-like charge carriers which behave like 
weakly interacting spinless fermions in the low-density limit L — N <C L. In this section, we 
use a reformulation of the Bethe ansatz equations due to Woynarovich || to construct an 
explicit mapping for arbitrary U and £5 <C 1 of the low- lying charge degrees of freedom of the 
ID repulsive-interaction Hubbard model onto spinless fermions whose mutual interactions 
vanish as £5 — > 0, and which obey a quadratic energy-momentum dispersion relation with a 
[/-dependent effective mass. We use this mapping to obtain asymptotically exact expressions 
for the charge stiffness and low temperature thermopower near the metal-insulator transition. 
Our calculation of the thermopower extends and makes more rigorous previous work of Schulz 
T0|1 which was based on a weak-coupling approximation, and implies that the transport is 



hole-like for < 1 - n < 

The energy eigenstates of a U > Hubbard chain of length L containing L — H electrons 
can be obtained by solving Eqs. ( p. 10] ) and Q2.11|) , with an appropriate choice of charge 



quantum numbers {/„} and spin quantum numbers {</<*}. The charge quantum numbers of 
the ground state are L — H consecutive integers (or half- integers), centered about the origin, 
and the spin quantum numbers of the ground state are (L — H)/2 consecutive integers (or 
half-integers), also centered about the origin (as discussed in Sec. 0). We wish to consider 
the low-lying charge excitations of the system, so we restrict our attention to those states 
whose spin quantum numbers {J a] are those of the ground state, and in which all of the 
pseudomomenta are real (states with complex pseudomomenta are separated from these 
by an energy gap ||). A hole Ih in the distribution of the charge quantum numbers (i.e., 
I n+ i — I n = 2, I h = I n + 1) corresponds to a charge excitation which, following Anderson p6| . 
we refer to as a holon. We wish to consider the case H <C L; since the I n are only defined 
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modulo(L), we can characterize the low- lying charge excitations of the system by the H holes 
{Ih, h = 1, . . . , H} in the distribution of the charge quantum numbers, rather than working 
with the larger set {/«}• Formally, we shall consider these holes in the distribution of the 
charge quantum numbers to exist even in the ground state of the system, where they will 
be consecutive integers (or half- integers) centered about L/2 (for convenience, we restrict 
Ih to the region 1 < Ih < L). We denote the corresponding holes in the distribution of the 
pseudomomenta by the set {kh}. Woynarovich M has derived a reduced form of the Bethe 
ansatz equations, valid in the large-L limit, which deals only with the parameters {Ih} and 
{kh} of these excitations, rather than the parameters of all the electrons. We extend the 
results of Ref. |J to include a magnetic flux $j = $^ = $ c through the ring, and point out 
that the reduced set of Bethe ansatz equations derived in Ref. |§ can be interpreted as an 
asymptotic Bethe ansatz for holons, which consequently is expected to become exact in the 
limit L — > oo for arbitrary H. We use this asymptotic Bethe ansatz to investigate in detail 
the limit H, L —> oo, with 5 = H/L small but finite. 

In the large-L limit, the energy of such a state is given in terms of the set {k h } by || 

H 

E(L -H) = E (L) - J2 £ M, (6.1) 

h=l 

where Eq(L) is the ground state energy at half filling ||, 

E (L) = -4 L r M»)Mu)du> (62) 

0V ; Jo uj[l + exp(uU/2)Y V ; 

and § 

eM--2cosfc-4 f°° J ^)cos(uj sin k)duj 

e c {k)- 2 cos k Aj Q u[1 + exp(fi)U/2)] ■ ^ 
The momentum of the state, defined modulo(27r), is |J 

H 

P = $ c - Pcih) + vr(L - N/2 + 1), (6.4) 

h=l 



where 



P ,(t)=H2r J -^ i "" . (6.5) 
Jo + exp(uU/2)\ 



Following Anderson ||46|| , we interpret —e c (kh) and — p c (kh) as the energy and momentum of 
a holon. The holons can not be regarded as noninteracting quasi-particles, however, since 
the kh are not free parameters, but are related to the set {Ih} by the equations |1 



H 



Lp c (k h ) = 2nl h + $ c + E (^> fcfc'), ( 6 - 6 ) 



h'=l 

where 
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ft n. k >) _ 9 /°° sin[a;(sinfc - sinfc')] cL; 

U ^ k )- 2 J ^[l + exp(^/2)] ■ (b ' 7j 



We point out that the set of equations ( |6.6D is formally equivalent to an asymptotic Bethe 
ansatz for the holons, G(fc, fc') being the effective holon-holon scattering phase shift [[|7j]. An 



asymptotic Bethe ansatz |$!| can be used to obtain the asymptotic form of the many-body 
wavefunction of a ID system with non- diffract ive scattering in the limit where the particles 
are widely separated — even if the exact wavefunction of the system is not of Bethe ansatz 
form. Moreover, the asymptotic form of the wavefunction is sufficient to determine the 
energy eigenvalues of the system exactly in the limit L — > oo, even at finite particle density 



35| . Consequently, we expect Eqs. (|6.1|), ( |6.4| ), and ( |6.6| ) to become exact in the limit 
L — > oo for arbitrary H. We emphasize that this conclusion follows a posteriori from the 
form of the equations, and is by no means obvious from the method || of their derivation. 
As a test of this conjecture, we use Eqs. ( |6.1|) and ( |6.6|) to calculate the energy of the empty 
lattice, which in this formalism is treated as the ground state of a system of L holons. In the 
limit H = L — > oo, Eq. ( |6T6| ) implies that the ho Ion pseudomomenta {k h } are equally spaced 
with density L/2ir, so that the sum over kh in Eq. ( |6.1|) can be replaced by L J^ n dk/2n. 
The resulting integral is equal to E (L), yielding E = 0, as expected for an empty lattice. 

Equations (|6.3Q and (|B~5| ) implicitly define an energy band e c {k{p)) for charge excitations, 
k(p) being the inverse of the function p c (k) defined in Eq. ( |6 . 5| ) . The holons are just holes in 
this energy band, which is full at N = L. Since de c {k{p)) / dp\ p=n = d 3 e c (k(p)) / dp 3 \ p=7T = 0, 
the holon energy near the zone boundary p = k = n can be written 

eMp)) ^r^r + „ + • • • , 6.8 

2\m*\ 4! 

where /i_ = e c (ir) is the chemical potential in the limit n — » 1~ 0, 



1 _ d 2 e c (k(p)) 



and 



A = 



m* dp 2 
d 4 e c (k(p)) 



4 ' (7r) <0, (6.9) 



p=7T 



^(tt)P 



dp 4 



p=TT 



[p'c( 7i )Y \p'A*)Y 



(6.10) 



Eq. ( |6.9| ) defines the effective mass m* of the holons near the zone boundary p8| , ^9| , |50 
which is negative, as is appropriate for a hole-like charge carrier. Explicit expressions for 
p' c (-K), s'c(tt), pi 3) (-K), and 4 4) (vr) are given in Eqs. (|B9])-( |BT2|) in Appendix B. 

Eq. ( p.6|) implies that the holon momenta differ from those of noninteracting spinless 
fermions by a term which vanishes as 5 — > 0; we write pt = p c {kh) = 2TrIh/L + Q c /L + Sph, 
where 5ph = L~ l J2h'=i ©(^/n kh>). Eq. (|6.7|) implies that 

4 In 2 

lim Q(k h , k h ,) = - [Pc(k h ) - p c {k h >)\ . (6.11) 

kh,k h i-^-K Up c y7T) 

For small 5 = H/L and low temperatures [ksT Ep ~ (7r5) 2 /2|m*|], we can use this 
expression for to calculate Sph' since the thermal average of 2irIh'/L is tt, the thermal 
average of 5pt (for fixed Ih) is 
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Inserting this result into Eq. ( |6.8|) , one finds that the shift of the holon momentum —5ph 
caused by holon-holon scattering leads to a fractional shift in the holon energy, measured 
relative to the chemical potential at half filling, given by 

-Se B 81n25 f -4 In 2 £8h U^O /fi10 . 

~~ ' (6.13) 



fi--e c Up' c (n) \-S\n26/U U -> oo. 

The effective holon-holon interaction, which should be considered as a correction to the 
hard-core repulsion of free spinless fermions, is thus attractive, and is negligible for £5^1 
(or U 3> t). While the explicit form of Eq. ( |6.13|) is only valid for £5 <C 1 and small 
excitation energies, one can show quite generally that \5ph\ < 8ln2 5/U, so that holon-holon 
interactions are also negligible at high temperatures when 5/U <C 1. 

The ground state of the system is obtained by choosing the set {Ih} to be consecutive 
integers (or half- integers) centered about L/2, i.e., by placing the ho Ions near the energy 
minimum at p = tt (with our convention for Ih, ph is restricted to the interval < ph < 2ir). 
Since Sph/ph ~ C(£<5) as 5 — > 0, the holon momenta are approximately equally spaced in 
the ground state when £5 -C 1. We can thus use Eqs. fl6.1|) , ( |6.6| ), and (|6.8|) to obtain an 
analytic expression for the charge stiffness near half filling [|5T|j, 



5 n 2 \6 3 2vr 2 ln2^ / j/ c (n) +p®(n) \ 

To 0(<5 3 ), the charge stiffness is equivalent to that of noninteracting spinless fermions in the 
energy band Q6.8Q , reflecting the approximate Galilean invariance which holds at low holon 
densities. The term proportional to 5 4 comes from holon-holon interactions. The term linear 
in 5 in Eq. (|6.14j) may also be obtained by a more conventional technique, which we describe 
in Appendix B. 

It is instructive to consider the weak- and strong-coupling limits of Eq. (|6.14j ). Making 
an asymptotic expansion of the ^/-dependent coefficients in Eq. QS.14p about U = 0, and 
using Eqs. (|3.7|) and ( |3.9|) , we obtain 



KmD c = ^ I 1 - \{^f + ^(^) 3 + ^(^) 4 ) > (6-15) 




which is equivalent to Eq. ( (4.5[ ), apart from the nonsingular multiplicative factor c/2t = 
1 + U/Airt. Note that the weak-coupling limit of the holon effective mass can be written as 
limjy^o \m*\ = A/c 2 , which is just the rest mass of the soliton in the equivalent sine-Gordon 
model. In the limit U — > oo, Eq. fl6.14Q becomes 



Jim D c = - (n5 - ^2 + 0(5 5 )) , (6.16) 

where we have made the hopping matrix element t explicit. This is just the Taylor series for 
(t/ir) shi7r5, the charge stiffness of noninteracting spinless fermions with dispersion e c (p) = 
—2tcosp, in accord with the well known behavior of the ho Ions in the limit U ^> t P,|37||. 
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The low-lying charge excitation eigenstates of a Hubbard ring with L+H' electrons follow 
from the above results and the particle- hole symmetry present at half filling || . The energy 
and momentum of such a state are given by Eqs. ( |6.1| ) and ( |6.4| ) with —e c {k h ) — > U — e c (kh), 
and the kh are determined by Eq. ( |6.6[) with $ c — > — $ c . For H' <C L, and for small 
excitation energies, these states can be described as the states of H' noninteracting spinless 
fermions (antiholons) near the bottom of the parabolic energy band, 

. C UHB (P) = U- eMv)) = U-^ + ^£ - + • • ■ . (6.17) 

2 1 m* I 4! 

The antiholons are thus characterized by a positive effective mass near p = ir. Note that the 
gap between the maximum of the lower Hubbard band Q6.8Q and the minimum of the upper 
Hubbard band fl6.17|) is U — 2//_ = 2 A, the Lieb-Wu charge gap [§||5||. The result ( J6.14p also 
holds for n — 1 + 5. We can thus write 

dn 2m* ' n ^ ' (6 8) 

where it is understood that m* > (< 0) for n > 1 (< 1). This expression emphasizes that 
the conductivity is "hole-like" when m* < 0, in that the Drude weight decreases when an 



electron is added to the system, and conversely, that it is "electron-like" when m* > |52 

The physical significance of the effective mass m* is revealed in the dynamics of a wave 
packet. The charge velocity v c is defined as the the group velocity of a wave packet composed 
of low energy charge excitations. For 5 Cl, we obtain 



de c {k{p)) 



dp 



P=tt(1—8) 



7l5 

r^i- 6 - 19 



If we apply an electric field E to the system, the wave packet will be accelerated []53 

^ = (¥) = - (- eE ) • ( 6 - 2 °) 
dt dp 2 \dt ) m* y 1 y 1 

The holons, having m* < 0, will be accelerated parallel to an applied electric field, and are 
thus hole-like, while antiholons are electron-like |54j . 



The transport coefficients which are conventionally used to determine the sign of the 
carriers — the thermopower and the Hall coefficient — couple to electrons, not holons, however, 
so the above picture is not the whole story. In ID, only the thermopower (Seebeck coefficient) 
S is available, which is defined in terms of the open circuit electric field E produced by a 
temperature gradient VT across the sample: 

E = SVT. (6.21) 



S is the entropy carried per unit charge by an electric current [ EE\ , and is ordinarily negative 
for electron conduction and positive for holes. In general, both the charge and spin entropies 
will contribute to S. However, the holon density of states diverges (v c — » 0) as 5 — » 0. Close 
to the metal-insulator transition, therefore, the entropy of the holons will be much greater 



than the spin entropy at low temperatures, and will dominate the thermopower |10,56 
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Furthermore, in this limit the holons can be treated as noninter acting spinless fermions in 
the energy bands ( |6.8|) and ( |6.17 ). So far, we have proven this only for the case where the 
spin wave function is in its ground state and there are no charge excitations into the upper 
Hubbard band [i.e., no excitations from the band (|6.8| ) into the band ( |6.17|) 1. Such charge 



excitations correspond to umklapp processes, and lead to a finite conductivity for T > |57 



however, their effect on the thermopower is negligible when fc#T <C 2A. The interaction of 
the holons with a thermal population of spin excitations can be described by adding a term 
L5p s h (T) to the right hand side of Eq. (|6.6| ); in Ref. H56fl , we show that for small 5 and low 
temperatures 5p s h (T) ~ 0(T5), so the holons are still effectively noninteracting when 5 and 
T are small. A standard formula for the low temperature thermopower |58|] can therefore 
be used to obtain 



S 



k 2 B Tm* 
IMl 2 



(6.22) 



which is valid for £5 -C 1 and ksT -C Ep — (tc5) 2 /2\m*\ in the absence of impurity scattering. 
The low temperature thermopower thus becomes large and positive as the metal-insulator 
transition is approached from n < 1, and has the opposite sign for n > 1. The small-f/ 
behavior of Eq. (|6.22j ) is in agreement with an earlier result JlO| obtained using a weak- 
coupling approximation. The entropy carried by spin excitations is important at higher 



temperatures and/or dopings, and is discussed in Ref. [56 . 



VII. CONCLUSIONS 

We have obtained an exact result for the charge stiffness of a Hubbard ring with U > 
and n — 1 in the large-circumference limit, which defines and yields an analytic expression 
for the correlation length £(£/) in the Mott insulating phase of the ID Hubbard model. We 
have shown that this correlation length also governs the exponential decay of the equal-time 
single-particle Green's function at n — 1 in both the weak- and strong-coupling limits, as 
well as the pairing correlations in the ground state of the attractive ID Hubbard model 
at n = 1. We remark that the strong-coupling expansion [Eq. ( p. lip ] which gave D C (L), 



G aa (L) ~ exp(— L/£), with = ln(U/t) + constant, must hold for the half-filled Hubbard 
model in the large-?/ limit in any spatial dimension. 

In the vicinity of the zero temperature critical point U — 0, n — 1, we have shown 
that the doping and system-size dependence of the charge stiffness scale with the correlation 
length £. The scaling function for the charge stiffness is expected to be universal ||, and 
thus should be applicable to other ID metal-insulator transitions. In addition, this scaling 
function appears to describe a certain class of ID magnetic phase transitions which do not 
involve broken symmetry, and which are characterized by a transition from a state with 
algebraic correlations to a state with a gap. In particular, Eqs. ( |2.3| ) and ( |4.2| ) imply that 
the spin stiffness D s = (L/2)d 2 E /d& 2 \q, s=0 (where $ s = ($ T - $|)/2) of the ID attractive- 
interaction Hubbard model with n = 1 and magnetization a = 2S Z /L has the form 

D a (S t ,L,-U)=Y±(fr,Z/L) (7.1) 

in the vicinity of the magnetic critical point U — 0, a — 0. The spin stiffness is related to 
the magnetic susceptibility \ by D s = (1/2tt 2 )x~ 1 (see Ref. [|]]). There is a similar mapping 
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between the metal-insulator transition occuring at n = 1/2 for ID lattice spinless fermions 
with nearest neighbor repulsion and the magnetic transition occuring at the isotropic point 
in the antiferromagnetic Heisenberg-Ising spin chain. Because of the universality of the 
scaling function for the metal-insulator transition, this magnetic transition is also expected 
to be described by the scaling function Y. 

One important consequence of scaling is that the finite-size corrections to D c are en- 
hanced near the critical point. While the scaling form ( |5.1| ) for the finite-size corrections 
can only be expected to hold for ID systems, it is plausible that the finite-size corrections 
to D c may be enhanced near half filling for smaller values of U in higher dimensional sys- 
tems as well, suggesting that caution is required in interpreting numerical calculations of 
o"(u;) on small clusters in this regime. Enhanced finite-size effects of this type could explain 
anamolous negative values for the Drude weight of the 2D Hubbard model on a 4x4 lattice 
with 14 electrons when U = 4 and 8 [59|. We do not expect similar enhancements of the 
finite-size corrections in the t — J model, since the constraint of no double occupancy fixes 

e = o. 

In the metallic phase of the model, the physical significance of the correlation length is 
that it defines the characteristic size of the charge- carrying solitons, or holons. We have 
shown that a reformulation || of the Bethe ansatz equations of the ID Hubbard model in 
terms of the parameters of the charge excitations only is formally equivalent to an asymptotic 
Bethe ansatz for holons, and have used this asymptotic Bethe ansatz to show explicitly for 
arbitrary U that the holons are equivalent to spinless fermions with mutual interactions that 
vanish in the low-density limit £5^0. We have used this mapping to obtain an expression 
for the low temperature thermopower near the metal-insulator transition, which implies 
hole-like transport for < 1 — n <C 

It has been argued in a series of papers by Anderson []6D| that the physics of doped 
Mott insulators is similar in one and two dimensions. It is therefore interesting to compare 
the metal-insulator transition in the ID Hubbard model with that observed in the cuprate 
materials in which high temperature superconductivity occurs, which are widely regarded to 
be quasi-two-dimensional doped Mott insulators ||. The thermopower of the doped cuprates 
is generically positive for hole doping 



3| and negative for electron doping 



35.66 



with a magnitude which increases drastically as the nominal concentration of doped carriers 
goes to zero, in qualitative agreement with Eq. ( |6.22| ). However, linear-T thermopower has 
not been clearly identified in metallic samples |)7| • F° r larger dopings, the thermopower has 



an unusual temperature dependence [ ]6T| , |6^|6^|6^|55] , |56] , |55| , |BU[ , and will be discussed in more 
detail elsewhere |56 . 

The optical conductivity near the metal-insulator transition in the doped cuprates 
70| , [7l|j72]| also shows qualitative similarities to that near the Mott-Hubbard metal-insulator 
transition in ID In particular, a linear growth of the Drude weight with doping 

is observed in the cuprates for small dopings |70| , |7T| , |72| , while the total integrated spectral 
weight in the optical conductivity up to ~ 3eV, which includes both the low-frequency re- 
sponse and the charge-transfer band (or upper Hubbard band) is roughly constant from 
the insulating antiferromagnet all the way through the superconducting phase [7(|[Fl]. The 
primary effect of doping in these compounds is thus to transfer spectral weight from high 
to low frequencies. This transfer of spectral weight is quite rapid, being essentially com- 
plete for a doping of ~ 0.25 [ fT0| , [7T| . Similarly, near the critical point of the Mott-Hubbard 
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metal-insulator transition in ID, we find that the total optical spectral weight 7rN to t is ap- 
proximately independent of doping, while the Drude weight tcD c grows linearly for small 
dopings, the transfer of spectral weight from high to low frequencies being essentially com- 
plete when £8 ~ 1. Of course, the fraction of the low- frequency spectral weight which is 
collapsed into the Drude peak at uj = is much greater in our ID model for several reasons, 
among which are the absence of an impurity potential, and the kinematic constraints which 
limit carrier-carrier scattering in ID. 
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APPENDIX A: CALCULATION OF f (J7) 



Here we derive the result ( |3.1| ) for the asymptotic form of D C (L) as L — > oo at n = 1 
using a technique JFjJ previously employed for computing finite-size corrections to energy 
eigenvalues in Bethe ansatz solvable models. 

We use Eq. (|2.4j) to calculate the charge stiffness from the dependence of the ground 
state energy on $ c . E ($ c ) may be obtained by solving the Bethe ansatz equations (|2.10|) 
and ( P- 11| ), and using Eq. ( |2.8|) . We consider only even N (= L), so that the ground state 
is a singlet (M = N/2), and nondegenerate. The charge quantum numbers are then: 



L 

-2 +1 



L 



L 

s 

2 



(Al) 



where s = j(L mod 4), and the spin quantum numbers are: 



M - 1 M-3 



M — 1 



(A2) 



In order to facilitate a solution of Eqs. ( |2.10| ) and ( |2.11 ) in the large- L limit, let us define 
the following functions |73| : 



p L (k) = k 



iV2tan- 1 ^ 



sin k 



U/A 



(A3) 



/ a x _ 1 V- o . _ x A - sin k n 1 ^ -iA-A^ 



Z7/4 



U/2 



With these definitions, the Bethe ansatz equations ( |2.10 ) and ( p.llj ) become, for $j 



Vh{k r , 



2nl n , $ c 



L ' 



(A4) 
(A5) 
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2nJ a 
L 



We also introduce the functions [f73 



p L {k) 



1 dp L (k) 
2tt rfA; 



(A6) 



(A7) 



^>4^' < A8 » 

In the limit L ^ oo with Af/L = n kept finite, the pseudomomenta {k n } in the ground state 
are distributed continuously on the real axis between the pseudo-Fermi points Q- and Q + 
with density p{k) = lim pi(k), and the spin rapidities {A Q } are distributed continuously 

L— >oo 

between — oo and oo with density //(A) = lim /i^(A). Equations (|2.10| ) and ( |2.11|) then lead 

L — >oo 

to the following coupled integral equations for p(k) and p(A) |74j]: 



27ip(k) = l + cosA; f°° dA— 8 ^ (A ) — , (A9) 

HK ' J-oo U 2 + 16(A-sinfc) 2 v ; 

2vm(A) + / dA' — , r = / dk — -f^ — . (A10 

m ; J-oo C/ 2 + 4(A-A') 2 JQ- f/ 2 + 16 A -sin A; 2 ' v ; 



where Q- and are determined by the conditions 
i-Q+ rQ+ 

IQ 



/ dkp(k)=n, / dkkp(k) = n$ c /L, (All) 
and is normalized to 

/oo 
dAp(A)=n/2. (A12) 
-oo 

The ground state energy is then 

E = -2L [ Q+ dk cos kp{k). (A13) 



E is independent of $ c at n = 1 in the limit L — > oo, as may be seen from Eqs. (|A9|) - 
( |A13| ): at half filling, the pseudomomenta span the entire Brillouin zone (— 7r, 7r), so the 



shift in the pseudo-fermi points Q_ and Q + caused by the flux $ c merely takes the Brillouin 
zone into itself, so that p(k), p(A), and Eq are all independent of $ c at n = 1. To obtain 
an expression for the ground state energy of the finite system suitable for an asymptotic 
expansion about L = oo, we express the sum over pseudomomenta in Eq. ( |2.8p as an integral 
by means of Eqs. (|A1[), (|A5|), and the Poisson summation formula, with the result: 



fir 00 

E ($ C ) = -2L dkcoskp L (k) exp{im[Lp L (k) - $ c + 2ns]}. (A14) 

J — ir ™_ „, 
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In the large- L limit, the dominant $ c -dependence of E$ comes from the terms with m = 
0, ±1. Moreover, for large L, Ph{k) may be approximated by p c (k) = lim pi{k) in the terms 

L— >oo 

with m^O. p c (k) may be obtained by integrating the expression for p{k) in Ref. ||, and 
was given in Eq. ( |6.5| ). Thus, using Eq. fl2.4j ), the charge stiffness may be written in the 
large-L limit as D C (L) = + D^ 5k \ where 

J-7T <5$ 2 



(A15) 



D (Sk) = _ r dksink s in[Lp c (k) + 2ns}. 

7T J-7T 



(A16) 



E)(Sp) describes the flux dependence of the ground state energy due to the change in the 
distribution of pseudomomenta as a function of $ c , while D^ Sk ' describes that due to a 
uniform shift of the pseudomomenta, 5k n = <& C /L. Both vanish as L —>■ oo. 
The integral in Eq. ( |A16| ) is dominated by the saddle point 



k = 7T + zsinh -1 ([//4), 
and may be evaluated by the method of stationary phase, yielding 
nm = {-l)* +1 LW{U/4) 

(L even), where 



\np»(k )/2\^ 



■exp[-L/e(C0]; 



oo 



(A17) 



(A18) 



1/£(U) = -ip c (k ) 



4 f°° \n(y + V2/ 2 - 1) 



U 



dy 



cosh(27ry /U) 



(A19) 



We next consider Eq. ( |A15|) for Dg p \ To find 5 2 pL(k)/5&l\$ c=0 , we first introduce 
coupled integral equations for piik) and //l(A). Combining Eqs. ( |A3|) , (|A7|) , and ( |A9|) , we 
obtain 



/oo 
dA 
-oo 



8^ L (A)-MA)] 
""[/ 2 + 16(A-sinfc) 2 



+ cos k 



dA 



8U 



U 2 + 16(A-sin£;) 2 



M 



Y,6(A-Ap)-f, L (A) 



(3=1 



(A20) 



The Poisson summation formula and Eqs. (|A2|) and ( |A6| ) can be used to rewrite Eq. (|A20[) 

as 



2ii[p L (k) — p(k)} = cos/c 



dA 



/oo 
dA 
-oo 



U 2 + lQ(A-sinky 



8U[p L (A) - p(A)} 
U 2 + 16(A-sin£;) 2 



exp{im[Lz L (A) + 2irs']}, 



(A21) 



m=— oo 
m^O 
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where s' = |[(M— 1) mod 2]. After similar manipulations, Eqs. (|A4|) , (|A"8|), and ( |A10| ) lead 
to 



27r^(A)-/i(A)] 



J17[>l(*)-p(*)] 



dk 



U 2 + 16(A-sinfc) 2 



00 dh! 4C/[MA') - MA')] 



f/ 2 + 16(A-sinA;) 2 



£/ 2 + 4(A-A') 2 

00 

^ exp{zm[Lp L (A;) - $ c + 2vrs]} 



m=— 00 
m^O 



<l\' ^77r~~~~7~~~~~ - f] exp{zm[ J L^(A / )+2 7 rs / ]}. 



f/ 2 + 4(A - A') 2 



(A22) 



m=— 00 
m^O 



Differentiating Eqs. ( |A2f|) and (|A22j ) twice with respect to $ c , and keeping only the leading 
order terms, we obtain the following coupled integral equations for 5 2 pi,(/c)/<5$ 2 |$ c= o and 
5 2 /xl(A)/5$ 2 |$ c=0 , valid in the large- L limit: 



2tt 



S 2 p L (k) 



5$; 



cos k I dA 



* c =o 



8U 



5Vl(A) 



t/ 2 + 16(A-sinfc) 2 5$ c 2 



<I>c=0 



(A23) 



2tt 



5Vl(A) 



<5$ 2 



+ / dA' 



*c=0 



4C/ 



5 2 MA') 



f/ 2 + 4(A - A') 2 
8U 



* c =o 
S 2 PL (k) 



U 2 + 16(A-sinA;) 2 5$ 



<I>c=0 



U z + 16(A — sin A;) z 
These equations can be solved by Fourier transforms, with the result 



5Vl(A) 



*c=0 



7T /-C 

— cos[Lp c (A;) + 2tts] / 

7T Z7T J-< 



exp[ico>(A — sin k)] duj 
2 cosh(u;C//4) ; 



(A24) 



(A25) 



*c=0 



^p(k')cos[L Pc (k') + 2ns]^%^ 

-TV In ok 



(A26) 



where Q(k, k') was defined in Eq. (6/?). Inserting Eq. ( A26 ) into Eq. ( A15 ), and evaluating 
the integral over k' by the method of stationary phase in the large- L limit, we obtain 



r\(5p) = (-l)^V 2 [l + (4/^) 2 ]V 2 exp[-L/e(^)] 
' " W(ftb)/2|Va 

da; e _IE tanh(x) Ji(4x/U); L — > 00 



(A27) 



(L even). The sign of is opposite that of .D^). As noted in Ref. 0], is negli- 

gible compared to in the large-?/ limit. In the small-?/ limit they are comparable in 

magnitude, but \D^\ > \D^\ V U. Combining Eqs. ( gig) and QA27D yields Eq. fl-TID, 
with 
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D(U) 



{U/Af - [1 + (?7/4) 2 ] 1/2 /o 00 ^e-nanh(a;)J 1 (4x/f/) 
(f J™dxe- x J (4x/U){(U/4) 2 + [1 + (f//4) 2 ]xtanh(x)}) 1/2 



(A28) 



By applying arguments similar to those above directly to Eq. ([A14 ), one can readily 
derive Eq. (|3T 



APPENDIX B: CALCULATION OF D c NEAR HALF FILLING 

Here we use a technique due to Haldane [[fj| to calculate the charge stiffness near half 
filling, verifying the term linear in 5 in Eq. ( |6.14| ). Eliminating /i(A) from Eqs. ( |X9| ) and 
( |A10|) yields an integral equation for p{k) alone fl76fl : 



(Bl) 



where Q(k,k') was defined in Eq. flOp . We note that Eqs. flATTD , (|AT3|) , and (|BTD , which 
determine the dependence of the ground state energy on $ c in the limit L — > oo, are identical 
in form to those considered by Haldane in his treatment of the spinless ID quantum fluid 



75fl . The result for the charge stiffness obtained in Ref. |75] is therefore applicable here, and 



is 



e 2< V 
2vr 



:b2) 



where = 1 — uj(Q) —uo{—Q) is a positive number which parameterizes the electron-electron 
interactions and 



2vrp(g) 



dke'(k)r(k) 



(B3) 



is the charge velocity. Here e(k) = —It cos k, and u;(A;) and r(k) are solutions of the following 
integral equations: 



2itu{k) = e(k,Q) - f Q dk' 9Q ^ — ou(k') 



(B4) 



2nT(k) 



de(k,Q) 
dQ 



-o dk 1 v ; 



(B5) 



Using Eq. fl6.7p for Q(k, k') and the fact that Q = tt at half filling, we see that v c vanishes 
at half filling because s'{ti) = and r(— k) = r(k) when Q = n. Also, uj(tt) = uj(—tt) = 0, 
so = 1 at half filling. Thus, for a small doping 5 away from half filling 



-6 



dD r 



dn 



-6 



dD r 



n=l~ 



dQ 



dQ 
dn 



(B6) 



n=l" 



From Ref. fTS] dQ/dn = [1 - w(g) + u(-Q)]/2p(Q), so (9Q/9n)| n=1 _ = l/2p(ir). After 
some manipulation, Eq. (|B6|) becomes 
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n tS ( 1 r . u d 2 B(k,Q) 
D c ~ 7 — 1 H / sin k -^n — - 

2np 7T 2 2tt J —-it dQ 2 



dk I . (B7) 

Q = 7T 



[27rp(7r)] 2 I 27T 7-. <9Q 
Inserting the result for p(ff) from Ref. into Eq. ( [B~T| ) and using Eq. ( |6.7| ), we obtain [51 



2 



V ^ 1 + exp(uU/2t)J / V •><> l + exp(u;*7/2*)/ v ; 

Note that the above expression for .D c may be rewritten using Eqs. (|6.9|) , ( p9|) , and (B10|) as 
D c ~ 5/ 1 2m* |, thus confirming the term linear in S in Eq. ( |6.14 ). Note also that v c ~ 7r<5/|m*| 
for 5 <C 1, in agreement with Eq. ( |6.19| ). 

For completeness, we give below the explicit expressions for p' c (7c), £"{Tf), ^■'(vr), and 



PM = 1-2I U ^ dUJ , , (B9) 
; Jo l + exp(uU/2ty y 1 



e»M = -2t(l-2r , ^toX (BIO) 



p w w=2 r ^ YTrr , (Bn) 



o 1 + exp(u;£V2t) 

+ exp(uU/2t) 



(B12) 



Using the Laplace transforms of Jo and Ji, and the Sommerfeld- Watson transformation, 
Eq. (p8|) may be rewritten in a form more suitable for an asymptotic expansion about U = 
as follows: 



D.*ts(* r * fa2 - 1)1/2[2c ° th2( ^ /lf> - 11 ) / ( r * (y2 r ir !^ 2 ■ Pis) 

I 17 A sinh(27rty/f/) / / \Ji y sinh(27rty/f/) / V ; 
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FIGURES 



FIG. 1. Ground state energy of a half-filled Hubbard ring versus magnetic flux {Tic/e)Q c from 
$ c = to 2vr. Solid curves: £ ($c) + N\U\/2 for U = -4; dashed curves: E Q ($ C ) for U = 4. (a) 
AT = L = 6 PS 1.5£; (b) N = L = 28 rs 7£, where £(4) ~ 4.06 is the correlation length at half filling 
for \U\ =4. 

FIG. 2. Plots of 7rD c versus f/L for (a) 5 = 0, (b) 5 = l/L, (c) 5 = 2/L, and (d) 5 = 4/L for 
systems with N = 60 (diamonds), = 80 (crosses), and ./V = 100 (squares), and 0.5 < U < 2, 
illustrating the scaling law, Eq. ( |4.2[ ). Note the different vertical scales, (e) For comparison, (b) is 
replotted with £ as the abscissa. 

FIG. 3. Total spectral weight 7riVtot (squares) and Drude weight 7T-D C (triangles) as a function 
of system size L for C7 = 6 and 5 = 0.2. The solid curves indicate L~ 2 behavior extrapolated to 
small L. 

FIG. 4. Plots of (a) A(5,U) and (b) B(6,U) as a function of f7 for <J = 0.05, 
0.1, 0.2, and 0.4 (top to bottom), where D C (L) = L> c (oo)[l + A{5,U)/L 2 + • • •] and 
N tot (L) = iV t ot(oo)[l + B(5, U)/L 2 H ]. Note the different vertical scales. 

FIG. 5. A plot of the function /+ = [A(5, U) -ir 2 /6}/£{U) 2 versus £(U)6. Solid curve: = 41; 
1.8 < U < 4.7. Triangles: U = 3; 5' 1 = 10, 15, 20, 30, 40, 60. Dotted curve: f+(x) = vr 2 /2 - vrV , 
the analytic result in the small-x limit. 
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